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We present measurements of radio emission from cosmic ray air showers that took place during 
thunderstorms. The intensity and polarization patterns of these air showers are radically different 
from those measured during fair-weather conditions. With the use of a simple two-layer model for 
the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation 
codes. This in turn provides a novel way to study atmospheric electric fields. 
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One of the important open questions in atmospheric 
physics concerns the physical mechanism that initiates 
lightning in thunderclouds [T]. Crucial to the answer 
is knowledge of atmospheric electric fields. Existing in 
situ measurements, from balloons or airplanes, are lim¬ 
ited due to the violent nature of thunderstorms. Fur¬ 
thermore, they are limited to balloon trajectories or per¬ 
turbed by the presence of the aircraft. Here we present a 
new method to probe atmospheric electric fields through 
their influence on the pattern of polarized radio emission 
emitted by cosmic-ray-induced extensive air showers. 

The main mechanism for driving radio-wave emission 
from air showers is that the relativistic electrons and 
positrons in the electromagnetic part of the shower are 
accelerated in opposite directions by the Lorentz force ex¬ 
erted by Earth’s magnetic field. This produces a short, 
nanosecond time scale, coherent pulse of radio emission 
mostly at megahertz frequencies. The emission generated 
by this geomagnetic mechanism is unidirectionally polar¬ 
ized in the Bvxb direction. Here, v is the propagation 
velocity vector of the shower and B represents Earth’s 
magnetic held US]. 

A secondary emission mechanism, contributing be¬ 
tween ~ 3—20% to the signal amplitude depending on 
distance to the shower axis and the arrival direction of 
the shower [SIS], results from a negative charge excess in 
the shower front. This consists of electrons knocked out 
of air molecules by the air shower. This also produces a 
short radio pulse but now polarized radially with respect 
to the shower symmetry axis. 

The emission from both processes is strongly beamed 
in the forward direction, due to the relativistic velocities 
of the particles. Additionally, the nonunity refractive in¬ 
dex of the air causes relativistic time-compression effects 


leading to enhanced emission from parts of the shower 
seen at the Cherenkov angle ] 710 . Interference between 
the differently polarized emission from both components 
leads to a complex and highly asymmetric intensity pat¬ 
tern ]5]. In contrast, time-compression effects do not 
alter the direction of the polarization vector of the emis¬ 
sion. The polarization pattern of the radio emission thus 
points predominantly in the Ovxb direction with a minor 
radial deviation. Strong atmospheric electric fields will 
influence the motions of the electrons and positrons in 
air showers. This is expected to be visible in the polar¬ 
ization patterns of the recorded emission [TO] . Therefore 
we analyze air showers recorded during thunderstorms. 

Data for this analysis were recorded with the low-band, 
10—90 MHz, dual-polarized crossed dipole antennas lo¬ 
cated in the inner, ^ 2 km radius, core of the LOw- 
Frequency ARray (LOFAR) radio telescope ]TT]. These 
antennas are grouped into circular stations that act as 
dishes for standard interferometric astronomical obser¬ 
vations. For the purpose of air shower measurements, all 
antennas are equipped with ring buffers that can store 
up to 5 s of raw voltage data sampled every 5 ns. A dedi¬ 
cated scintillator array, LOfar Radboud air shower Array 
(LORA), is located at the center of LOFAR to provide 
an independent trigger whenever an air shower with an 
estimated primary energy of > 2 x 10^® eV is detected 
]T2]. When a trigger is received, 2 ms of raw voltage data 
around the trigger time are stored for every active an¬ 
tenna. 

These data are processed in an off-line analysis [TJ] 
from which a number of physical parameters are ex¬ 
tracted and stored. These include the estimated energy 
of the air shower (as reconstructed from the particle de¬ 
tector data), the arrival direction of the air shower (as 
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reconstructed from the arrival times of the radio pulses in 
all antennas), and for each antenna polarization informa¬ 
tion in the form of the Stokes parameters: I (intensity), 
Q, U and V. The orientation of the polarization vector is 
reconstructed from Stokes Q and U [5]. 

Over the period between June 2011 and September 
2014, LOFAR recorded a total of 762 air showers. The 
complex radio intensity pattern on the ground of almost 
all measured showers can be well reproduced by state- 
of-the-art air shower simulation codes M- These codes 
augment well-tested Monte Carlo air shower simulations 
with radio emission calculated from first principles at the 
microscopic level [El IE]. In this analysis, we use the 
CoREAS plugin of CORSIKA [T7j with QGSJETII [E] 
and FLUKA [Ej as the hadronic interaction models. It 
was previously found that the exact shape of the inten¬ 
sity pattern depends on the atmospheric depth where the 
number of shower particles is largest, Amax, and that the 
absolute field strength of the radio emission scales with 
the energy of the primary particle. 

The radio footprints of 58 of the 762 air showers are 
very different from those predicted by simulations. Of 
these, 27 air showers have a measured signal-to-noise ra¬ 
tio below 10 in amplitude — too low to get a reliable 
reconstruction. The polarization patterns of the other 31 
showers differ significantly from those of “normal” fair- 
weather air showers. This can be seen in the middle and 
bottom panels of Fig. where the polarization direction 
is clearly coherent (i.e., nonrandom) over all antennas but 
no longer in the expected Bvxb direction. In addition, for 
some of these showers the intensity of the radio signal at 
low 10—90 MHz frequencies is strongest on a ring around 
the shower axis with a radius of approximately 100 m 
(see also Fig. [^. This “ring structure” in the intensity 
pattern is not present in normal fair-weather air showers 
that all lack rotational symmetry in the intensity pattern 
and instead show a single maximum that is displaced in 
the BvxB direction from the shower axis [Ml[20]. Twenty 
of these 31 showers occur within 2h of lightning strikes 
recorded within ~ 150 km distance from LOFAR by the 
Royal Dutch Meteorological Institute. Given the similar¬ 
ity of the polarization patterns of the remaining showers 
where no lightning strikes were measured, it is plausi¬ 
ble that at these times the atmospheric electric field was 
also strong albeit not strong enough to initiate lightning. 
An electric field meter has since been installed at LO¬ 
FAR that will provide independent verification for future 
measurements. 

For the shower in the middle panel of Fig. recorded 
during thunderstorm conditions, the pattern is unidirec¬ 
tional for the entire footprint. A second more compli¬ 
cated type is depicted in the bottom panel. Here, the 
pattern is more “wavy”. The analysis presented here fo¬ 
cuses on an air shower of the first type where also a ring 
structure is visible and a strong signal is measured by 
the LORA particle detectors. All air showers of this type 
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FIG. 1. Polarization as measured with individual LOFAR 
antennas (arrows) in the shower plane for three measured air 
showers. LOFAR antennas are grouped into circular stations, 
of which seven are depicted. The expected polarization direc¬ 
tion for fair-weather air showers is indicated with “normal”. 
The position of the shower axis, orthogonal to the shower 
plane, is indicated by the intersection of the dashed lines. 
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can be reconstructed with similar accuracy. For showers 
of the wavy type a more complex analysis is currently 
being developed. 

We propose that the influence of atmospheric electric 
fields on air shower radio emission can be understood in 
the following way. 

The electric field, in the region of the cloud traversed 
by the air shower, can be decomposed into components 
perpendicular Ej^ and parallel Ey to the shower sym¬ 
metry axis. The perpendicular component of the field 
changes the net transverse force acting on the particles 

F = g(Ej.-Hv X B). (1) 

This changes both the magnitude and the polarization of 
the radiation that follow F. 

During shower development the air shower particles 
lose energy. The parallel component of the atmo¬ 
spheric electric field partially compensates this energy 
loss. Therefore, the total number of particles within a 
given energy range in the shower increases. Because the 
fractional gain of energy is greatest for lower energy par¬ 
ticles, these are the most affected. However, low-energy 
particles do not contribute much to the total radio emis¬ 
sion because they lag behind the shower front and their 
emission is not coherent for frequencies above 10 MHz. 
Thus, it is the perpendicular component of the electric 
field that determines the measured intensity and polar¬ 
ization direction. 

In order to test these hypotheses, atmospheric elec¬ 
tric fields were inserted into CoREAS air shower simula¬ 
tions. By the comparison of fields acting purely parallel 
and purely perpendicular to the shower axis it was found 
that the effect of Ex on the radio emission is indeed 
much stronger and will dominate in most shower config¬ 
urations where both components are present. This will 
be discussed in greater detail in a forthcoming publica¬ 
tion. 

Having understood the basic effects of atmospheric 
electric fields on air shower radio emission we proceed 
with a full reconstruction of LOFAR measurements. We 
follow the method developed by Buitink et al. [Uj to fit 
CoREAS simulations to LOFAR measurements. An at¬ 
mospheric electric field is inserted into the simulations 
with the perpendicular component chosen such that the 
net force is in the measured average polarization direc¬ 
tion (as indicated in the middle panel of Fig. [^. The 
parallel component is set to zero since its influence on 
the received radiation intensity and polarization pattern 
is negligible. 

The simplest electric field configuration that can re¬ 
produce the main features in both the measured inten¬ 
sity and polarization patterns is composed of two electric 
field layers. The upper layer, with strength lE^^I, starts 
at a height hjj above the ground and extends down to a 
height Hl bX, which the lower layer starts, the direction 
of the net force changes by 180° and the field strength 


decreases to |Ex|. Two layers are needed because with 
one layer the ring structure seen in the measurements is 
not reproducible. 

In Fig.j^the reconstruction is shown for the air shower 
for which the polarization pattern is depicted in the 
middle panel of Fig. The reconstruction is optimal 
for hjj — 8 km, = 2.9 km, |E( 7 | = 50kVm“^ and 
|Ex|/|Et 7 | = 0.53. For these values y^/ndf = 3.2 as 
obtained for a joined fit to both the radio and particle 
data. A perfect fit of y^/ndf « 1, as is often found 
for fair-weather showers, is likely not attainable with a 
simplified electric held model. However, all the main fea¬ 
tures of the intensity and polarization pattern (namely 
the overall polarization direction and ring structure) are 
already correctly reproduced. 

The fit quality is sensitive to changes in the relative 
field strength and hx as well as Amax- This can be seen 
in Fig. where each parameter is varied while keeping 
the others fixed at their optimum values. This fixing is 
not possible for Amax in the CORSIKA software, because 
it is a an outcome of the simulation rather than an input 
parameter. Therefore, simulations were selected where 
Ainax differs by no more than 20gcm“^. The fit qual¬ 
ity reaches its optimum value for hjj = 8 km and is not 
sensitive to a further increase. This is expected because 
above this altitude the air shower is not yet fully devel¬ 
oped and there are relatively few particles contributing 
to the emission. 

For fair-weather air showers the measured radio inten¬ 
sity is related to the simulated values through a con¬ 
stant scaling factor |14] given the energy of the primary 
particle. This energy is derived from the particle den¬ 
sity on the ground, as measured with LORA, combined 
with the information on Amax, as determined from the 
radio fit. For the air shower measured during thunder¬ 
storm conditions the measured intensity is higher than 
the normally expected value, as the absolute electric field 
strength influences the radio intensity. However, the sim¬ 
ulated intensity increases only until the atmospheric elec¬ 
tric field strength reaches |Ec/| > 50 kV m“^. When the 
field strength is increased further the radio intensity stays 
constant. This saturation of the radio intensity appears 
to be related to the coherent nature of the emission but 
is still under investigation. 

Measuring radio emission from extensive air showers 
during thunderstorm conditions thus provides a unique 
new tool to probe the atmospheric electric fields present 
in thunderclouds. Unlimited by violent wind conditions 
and sensitive to a large fraction of the cloud this tech¬ 
nique may help answer the long-standing question “how 
is lighting initiated in thunderclouds?” It has been sug¬ 
gested by Gurevich et al. [2T1[22] that cosmic-ray-induced 
air showers in combination with runaway breakdown may 
initiate lightning. If this is indeed true then LOFAR with 
its combination of particle detectors and radio antennas 
is well positioned to measure it. 
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FIG. 2. Radio intensity pattern during a thunderstorm. Top: 
the circles represent antenna positions. Their color reflects 
measured pulse power. The best-fitting CoREAS simulation 
is shown in color scale in the background. Where the colors of 
the circles match the background a good fit is achieved. Bot¬ 
tom: measured (circles) and simulated pulse power (squares) 
as a function of distance to the shower axis. 
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FIG. 3. Sensitivity of the fit quality to variations in the atmospheric depth of shower maximum Xmax (left panel), the relative 
field strength (middle panel) and the field reversal altitude h 2 (right panel). The optimal proton simulation is the same for all 
plots. The electric held strength, in the upper layer, for all simulations is lEy] = 50kVm“^. 
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